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'sport on the Progress made in the Calculations for aNeic ’Method 
of treating the Lunar Theory. By Sir Gr. B. Airy, K.C.B., 
Astronomer Royal. 

In the Monthly Notices of this Society for 1874, January, I 
explained the form which I proposed to give to a New Method of 
treating the Lunar Theory, and I stated that I had actually 
commenced some of the numerical developments of which the 
new method mainly consists. In my usual Annual Report to 
the Board of Visitors of the Royal Observatory, I mentioned (in 
the ordinary course of report) that I had arranged such a 
method of treating the Theory, that the principal part of the 
work was purely numerical, and that a computer had been 
engaged, at my expense, in proceeding with some parts of the 
work. The Board spontaneously passed a resolution to the effect 
that this undertaking might properly be considered as a part or 
extension of the astronomical calculations made at the Observa¬ 
tory, and that in their opinion it might properly be undertaken 
at the public expense. I submitted this opinion of the Board of 
Visitors to the Admiralty, through whom it was placed before 
the Treasury; and, after explanation of the nature and object of 
the proposed work, I received authority from the Treasury for 
incurring the expense required for carrying on the computations 
through the current Financial year. I have subsequently sub¬ 
mitted an estimate for progress of work during the year 1875- 
1876. Two computers have been employed, the aggregate of 
their times of employment to 31st December 1S75 being about 10 
months. 

We have had painful evidence of the difficulties that may 
arise in utilising a laborious investigation, when by some sad 
chance the possibility of its completion by the original proposer 
has been suddenly terminated, and scarcely sufficient explanation 
has been left to enable successors to enter into the remaining 
steps in the same spirit as that which guided the primary opera¬ 
tions. It has appeared to me, therefore, desirable in prosecuting 
a new theory (which, though far simpler than preceding theories 
in the intellectual labour demanded, is still sufficiently extensive 
and sufficiently novel in form to require some study before a fol¬ 
lowing astronomer could actively take it up), to place on record 
such a description of the method of proceeding, and, from time 
to time, such a statement of the advance of the work, as will 
enable a successor, if necessary, to carry on the operations 
without much difficulty to their close. 

1. Referring to my paper “ Proposed New Method &c.” 
already cited, paragraph marked VI., it will be seen that the 
first step is to convert Delaunay’s sexagesimal angles, expressed 
in seconds of arc, into equivalent expressions in decimal fractions 
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1 of radius; that these expressions are to he limited to the seventh 
j decimal place, of which the unit is 10-7; and that, this limit 
■being understood throughout, trouble and confusion will be 
^voided by omitting decimal cyphers preceding the effective 
[figures. To convert the seconds of arc into decimals of radius, 
[it is necessary to multiply the number of seconds by 
-0-000004848137. Therefore, multiplying by 100000, and con¬ 
sidering the final figure of Delaunay’s coefficients to be the 
second decimal of the sexagesimal second, we have to arrange a 
method, applicable to all coefficients, by which 10000000 maybe 
converted into 4848137. The method adopted is one analogous 
to that known in books of arithmetic by the term “ Practice 
it never leaves a doubt on the value of the leading figure, and it 
employs no unnecessary figures. By trial, the following steps 
are adapted for this case ; every succeeding number being formed 
from the preceding number by use of the fraction in the margin : 


10000000 


+ ~ = -t- 5000000 


IOO 



— 150000 

I 

90 



— 1667 

I 

8 ^ 



— 208 
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20 ~~ 
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IO 
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H- 
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+ 

5000012 

-15187s 


~ 15187s 

+ 4848137 


In the application to actual coefficients, the last term is never 
wanted, and the preceding term only once. Of the numbers, 
which I now proceed to exhibit as specimens, the first is the 
argest coefficient in longitude, the second is one taken from 
the middle of the series :— 
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In succeeding parts of the work, several thousands of multipli¬ 
cations are performed in the same manner; the numbers to be 
multiplied being arranged in series, of which every term has the 
same multiplier; and the marginal fractions for that multiplier 
being determined by trial. 

In some cases it is necessary to perform division ; thus, in 
order to reduce the expression for parallax from the form 

Parallax = 3422700 4- series of terms 

to the form 


— = {1 + new series of terms}, 

it is necessary to determine by trial a process which will convert 
34227000 into 10000000. The marginal fractions are found to be 

1 1 1 1 1 

3’ ~ 9 ’ ~io’ ~ 9 ? ~~ 3 ° e 

2. The number of coefficients of cosines of different argu¬ 
ments in the expression for Parallax, adopted from Delaunay, 
and converted by the process above described, is 100. The 
number of coefficients of sines of different arguments in the ex¬ 
pression for Longitude, is 170; and that in the expression for 
Latitude, is 142; which are similarly converted. By subsequent 
operations, the numbers of coefficients in squares, &c., of these 
series are raised to 148, 170, 180. The smallest value retained 
is 1 (in my notation), or a number not less than its half. 
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!;! 3. The great mass of operations consists in multiplying 

together two series, such as 

[co[ 

!s! A + B cos fp] + C cos [q| + D cos |B| + etc. 

[loJ 

[Sj a + b cos |p] + e cos [5] + d eos |"r] + etc. 

The coefficients of the upper line are all multiplied by a (by the 
process above described); then by b, and so on. The general 
product is 


A«+ B«cos[P] + C«cos[q] + etc. + A5cos [p] +^cos |P +p[ + ^cos |P—p| 


+ ~cos ]Q + pj + cos |Q —p| + etc. + Ac cos fg] + ^ cos fFT$"| 
+ ^cos |P— g| + ^cos |Q, + g| + ^ cos |Q — jI + etc. 


In the great majority of instances, the new arguments formed in 
these multiplications are to be found in each of the series which 
are multiplied together ; and the form of the series produced by 
the collection of these products is the same as that of the first 
series, with the addition, perhaps, of a few small terms. In the 
primary series for Parallax, Longitude, and Latitude, the co¬ 
efficients are arranged in the order of descending magnitude, and 
this order is little disturbed in the operations of multiplying series 
together. 


4. The coefficients of- —1 are represented by the ordinal 
numbers (written in red ink) 1, 2, 3, &c., to 100 ; those of 
e-) by 201, 202, &c., to 348, the small figures following cen¬ 
turies being attached to the same cosines, as far as they go: 
those of (H' by 401, 402, &c., to 535, the small figures fol¬ 
lowing centuries being attached to the same cosines, &c. The 
coefficients of latitude are represented by 2401, 2402, &c., those 
of longitude by 5001, 5002, &o.; and the coefficients of products 
and other functions formed from them, by different series of 
similar order, the first always being a century increased by 1. 

5. When powers of - — 1, to the fifth power, have been 


formed by successive multiplications by the series for 1, 

(forming a new series for each power, of which the terms are 
distinguished by new centuries of red-ink numbers), all the re¬ 
quired powers of - are formed by the binomial formula : 
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r + n I - 


+ etc. 


AVhen powers of Latitude have been formed by successive multi- 
Ijlication by the series for Latitude, the expression for sine 

3 Latitude and the expressions for the various powers of cosine 
"["Latitude i are formed by means of the appropriate formulas in 
powers of Latitude. 

6. For the formation of —, it is to be remarked that every 

dt 

argument is formed by the sum of multiples of four simple 
arguments, namely, 


D = Excess of Moon's mean longitude above Sun’s mean longitude. 
f — Moon's mean argument of latitude. 

I - Moon’s mean anomaly. 

L = Sun’s mean anomaly. 


The proportion of the daily increase of each of these arguments 
to the daily increase of the Moon’s mean longitude is taken as 
adopted or found by Delaunay, and the unit of time is supposed 
to be such that the Moon’s mean motion in longitude = 1. Thus, 
for the term 

sin j2D + 2/ + 2 ~Ii|, 

the sum of the proportions of daily increases of terms under the 
bracket to daily increase of Moon’s mean longitude is found to 
be 47751 ; therefore 


^ sin 1 2 D + 2/ + l —L [ 


= cos 1 2D + 2 /+f—Lj x 47751 x — (daily increase of Moon’s mean longitude) 

Qj o 

= 4775 1 x cos I2D + 2/+ Lj. 

The corresponding term in ~ is found by multiplying this quan- 

Oft 

tity by the coefficient of the corresponding term in v. Succeed¬ 
ing, differentiations are effected on the same principle. 

7. The methods explained to this point apply to the expansion 
of functions on the basis of elements assumed by Delaunay, as 
representing the Moon’s motion at a standard modern year (say 
1900) ; without any consideration of the corrections which are 
due to the coefficients in those expansions, in order to produce 
a correct theory, on the supposition that the elements of the 
Earth’s orbit round the Sun are stable; and without any considera¬ 
tion of the corrections which are due to the slow changes of the 
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Elements of the Earth’s orbit, which produce secular terms in the 
I^Jfoon’s motions. Notice of these corrections will be taken here- 
1 after. At present I may witb propriety record the progress made 
;|n the limited expansions to which I have first alluded. 

8 . The following quantities (each expanded into a series of 
jipumerous terms with numerical coefficients and periodical 
“trigonometrical functions) are computed, all in duplicate, and 
some of the earliest in triplicate. 


Subject of Expansion. 

Modified equivalent of 

a , a 
- and — — i 
r r 

(?-■)’ - 

c-y. 

c-y - ••• 

c-y ••• 

d-y . 

(if . 

c-y . 

r 

a 



Equivalent for Latitude 
(Latitude) 2 ... 
(Latitude) 5 ... 


Reddnk Nos. 

I— IOO 
201— 348 

4 °*— 53 1 
601— 725 
801— 826 
1001—1148 
1201—1348 
1401—1548 
1601—1748 

I— IOO 

1801—1948 
2001—2148 
2201—2348 

24OI—2542 
2601—2726 
2801—2959 


(Latitude) 4 ... ... ... ... ... 3001—3128 

(Latitude) 5 ... ... ... ... ... 3201—3246 
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Subject of Expansion. 

Ked-ink Nos. 

(Latitude) 6 ... 

3401—3434 

sin [X] 

3601—3754 

cos \T\ 

3801—3927 

(cos [Hr. 

4001—4127 

(cos rri) 3 

4201—4327 

(cos nn y 

4401—4527 

COS 1 • -- 

it 

4601—4790 

( tM T*;)* . 

4801—4995 

Equivalent for Longitude 

5001—5182 

dv 

dt 

5201—53S2 

IK- m-;)* ••• . 

5401—5645 

Cfh - 

5701—5945 

. 

6001—6191 

. 

6201—6450 

(Velocity) 2 , or sum of last quantity + last but two 

6501—6750 

(Velocity) 2 

6801—7050 


(The velocity here mentioned is the resolved velocity 
parallel to the plane of the Earth’s orbit.) 

The duplicate calculation of the quantity last recorded is not 
yet completed. 

The calculation of the following quantities is not yet 
begun 




T 

a 


• sin |T . 


d ? (r . — \ 

S V “°) 

(;)' »»m• 


X 2 
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■ :■ These complete the list of all the quantities not multiplied by 
jtjhe disturbing force which enter into the Equations of the Lunar 
theory. 

Some of the computed numbers, applying to equations of 
|!png period, may require further revision. 

|®| 9. I now proceed to remark on the state of the “ symbolical 

“corrections ” mentioned in article (2*) of my former paper, on 
whioh the corrections of Delaunay’s coefficients, as depending 
on the elements of the Earth’s orbit in a given year (say 1900), 
are based ; and from which the secular equations of the Moon’s 
motion, depending on the slow changes of those elements, are to 
be deduced. For the present, I limit my details to the symbolical 
corrections of the coefficients of the periodic terms. The nota¬ 
tion used to express these symbolical variations is, the red-ink 
number defining the coefficient (as above described) surrounded 
by a ring. This applies to every coefficient in the first instance, 
and to every subsequent coefficient, except where, by virtue of 
multiplication by small numbers, it is made too small to possess 

any value. Thus the series for - — 1 is expressed by 

(red-ink 2 + (2^)) • cos [ 7 ] + (red-ink 3 + (3)) .cos I2D — 1 \ + etc.; 


the series for the longitude-terms added to mean longitude is 
(red-ink 5002+ ( 5002 )).sin [ 7 ] + (red-ink 5003 + (5003) j.sin|2D —£| + ete.; 


the series for latitude is 

(red-ink 2401 + ( 2401 )). sin |/| + (red-ink 2402 + ( 2402 )). sin |/+ 1 \ + etc., 

where “red-ink 2,” “red-ink 3,” “red-ink 5002,” “red-ink 
2401,” &c., are numerical quantities, adopted or computed as 
above-mentioned; but 

(2), (3), (5002 ), ( 2401 ), etc. 

are still preserved in the symbolic form. 

10. As the further treatment of these symbolical terms tends 
to increase their number and the complexity of their combina¬ 
tions, it is necessary from the beginning to adopt some rule as 
to the anticipated magnitude of the corrections which Delaunay’s 
coefficients may require. Now I think I may safely assume 
that no correction greater than 2 //, oo, or 100 of my notation, can 
result from the investigation. In this case, it is useless to retain 
any product by a number smaller than o‘ox, or (with allowance 
for turning of the figure) smaller than 0^005, or smaller than 
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50000 of my notation. In the improbable case of an error ulti- 
j mately showing itself with magnitude greater than 2"‘oo, the 
; only effect of our limitation would be, that in the equations there 
; is occasionally a loss of some quantity amounting numerically to 
; 1 in the last figure. 

; 11. The application of these considerations will be best 

understood from an example. The first instance of multiplication 
is that of the formation of 


0 “ l ) multiplying 0 -1) by 0 -1); 
or 


(r.i. 2 + (?)).cos |: 4- (r.i. 3 + ( 7 )). cos 2D — 1 \ + (r.i. 4 +( 7 )). COS } 2 l>{ 

+ (r.i. 5 + @).cos ,> 27 ! + etc.; 

multiplied by 


(’•“+©) .cos [ 7 ] + (] r.i. 3 + ( 7 )) -COS J ^ j + (r.i. 4+ (4)).cos 2D| 

♦('**+©) . cos 12 1 j 4- etc.. 


which will produce 

(r.i. 201 + (201)^ .cos [ol 4* (rA. 202 + (202. cos m + (r.i. 203 + (203)). cos 12D- 7 ] 
+ (r.i. 204 + (204)) . cos ; 7 d| + (r.i. 205 + ( 205 )). cos [77] + etc. 


Now, cos [oj is produced by the combination of cos [ 7 ] with 
cos [ 7 ], of cos I2D — 1 \ with cos |2 D — 1 \, of cos [2 D| with cos 1 2D], 
of cos [ 77 ] with cos [ 77 ], etc., and its coefficient, or r.i. 201 +(201), 
will therefore be 

7 2) 2 + r.i. 2 x ( 7 ]), 

+ \ ( r.i. 3) 2 + r.i. 3 x 

+ \ (r.i. 4)' + r.i. 4 x (7) , 

+ “(™. 5) 2 + r.i. 5 x 
+ etc. 


Now r.i. 2 = +545094 ; r.i. 3 = +99812 ; r.i.4 ss - 82321; 
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! r. i. 5 = 4-29794. Here then it appears that r.i. 2, r.i. 3, r.i.'4, 

■ are above our limit; but r. i. 5 and all succeeding numbers 
;ij which are still smaller) are below it. The first three, then, of 
;|jhe products of symbolical variations just given, are the only 
ijlines to be retained for (201). Striking off the last five figures 

loo I 1 

^we have 

(201) ^ + 5 x @ + l x(s) + ixQ), 


where the unit of the factors is 


unity 
100 ' 


or io~ 2 . 


Now consider the factor of cos j /1 in the term 

( r . i . 202 + ( 202 )) .COS [7]- 

Here, cos | 7 ] may be produced by the combination of cos [ 7 ] 
with cos | 77 ]j of 00s | 2 D —Z| with cos | 2 fj], of cos j 2 D | with 

cos | 2 3 D — Z|» of cos | 2 I) — 1 with cos |' 2 D — 2V1, and of others 

which are at once seen to be insensible. Hence, the value of 
(r.i. 202 +( 202) ) is 

+ (r.i. 2 + ( 7 )) x (r.i. 5 + ( 7 )) > 

+ (r-i- 3 + (T)) * ( r - L 4 + (4)) » 

+ (r.i. 4 + ©) x (™* 6 + (*))» 

+ (r.i. 3 + (?)) x (r.i. I 9 + (^9)); 

and the value of (202) is f . 


+ r.i. 2 x 

+ r.i. 5 x @, 

+ r.i. 3 x ( 4 ^) 

+ r.i 4 x Q), 

H- r.i 4 x 

+ r.i. 6 x 0, 

+ r.i. 3 x (g) 

+ r.i. 19 x (3). 


The only factors which are sensible are r. i. 2, r. i. 3, r. i. 4 ; 
and therefore 
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o ^ 1 
001 
CM 1 

(202) = + 5 * © + I x (4) + 1 x © + 1 x © + 1 X (g), 


or iO‘ 


Inhere (as before) the unit of the factors is - Uni ^ , 

■ .100 

| When the series for — 1^ is multiplied by another series, 


a similar process must be used for expressing the symbolical 
variations of the coefficients then newly formed, in terms of 

(201) , (202) , etc. • and of the variations of terms of the other 
series; then the values found above are to be substituted for 
(201), (202 ), etc., and finally the symbolical variation of every 
coefficient, without exception, will be found in terms of 


©), ©• (T). etc . (2401), (2402), etc., and (5001), ( 5002 ), etc. 


12. The symbolical variations are formed as far as 

(cos |71 X y 

in the table of article 8, but may require further revision. 

13. To complete the exhibition of' symbolical variations, it is 
necessary to prepare the variations produced by the symbolical 
variations of the arguments. Thus, for the term 

r.i. 8 x cos ] 2 D-/-L|. 
the complete variation is 

(?) . cos | 2 D —Z —L | - r.i. 8 . sin | 2 D-Z-L| x ( 2D-I-L ) . 

In the ordinary process of multiplying series by series, this opera¬ 
tion is remarkably simple; it is not necessary to go step by step, 
as in the preparation of the variations of the coefficients in 
article 11 ; but the variation depending on the argument will be 
found correctly from the argument standing at the end of the 
operations (as will easily be seen by trying a step-by-step opera¬ 
tion). A trifling complexity is introduced in these terms which 
bear the sign of differentiation, because the symbolical variations 
themselves are subject to differentiation. 

14. As the variation of arguments produces its effect in the 
grand final equation only where it is multiplied by a coefficient 
which is usually small, and as it is particularly in these variations 
that we must look for the terms which rise to importance in the 
shape of secular terms, we must be prepared to accept much 
larger values for the magnitude of the terras to be produced by 
these symbolical variations of arguments. I have not yet decided 
on the limit to be adopted. 
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!;! 15. No step of calculation, has yet been made in reference to 

■Sihe symbolical variation of arguments. 

|m| 16. No calculation is yet made for that side of the final 

■|equation which contains the disturbing force as a factor. 

lo 

T 1 

ly! In any future Report which it may appear desirable to place 
before the Society, there will be no occasion to repeat any part 
of the explanation of methods given here, and it may be antici¬ 
pated that the statement of the step which the operations may 
have reached will be very brief. 

Royal Observatory , Greenwich, 

1875, February 6. 


Observations of the Satellites of the Planet Uranus, made at Birr 
Castle , during the years 1872, 1873, and 1874. 

By the Earl of Rosse. 

The following observations of the Satellites of the planet 
Uranus , made with the reflector of six-feet aperture at Birr 
Castle, though for several reasons not as satisfactory as might 
be desired and might be obtained with that instrument, may 
perhaps be of some value as a contribution to knowledge in a 
department of astronomy in which so few have as yet worked. 

The shortness of the time (40 minutes) available each night 
for the observation of the planet with the 6-foot instrument, the 
atmospheric disturbance so often a source of annoyance in using 
so large an aperture, the unsuitability of the micrometer for the 
work, and to a certain extent also the almost unavoidable deteri¬ 
oration of the polish of the mirrors in the course of working, 
which probably increased the diffused light round the planet, 
combined to detract from the accuracy of the observations, and 
to render the detection of the inner satellites difficult. 

My late assistant, Dr. Ralph Copeland, has compared the 
observed places of Oberon and Titania with those which he has 
calculated from Dr. von Asten’s definite circular Elements given 
in Yol. xviii., No. 5, page 26, of the Memoirs of the Imperial 
Academy of Sciences of St. Petersburg. The differences between 
calculation and observation are given in the column headed 
(C. — 0 .). The observations made by Dr. Copeland and myself 
are distinguished respectively by the letter C. or R. 
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